* 


Semi-Annual Status Report 



($s() /J/!^ 

cjy A? / ^ 

/ V - 55 '- <?<- 




"The Absolute Radiometric Calibration of the 
Advanced Very High Resolution Radiometer" 


97/f> & 


P. N. Slater 
P. M. Teillet 
Y. Mao 


Optical Sciences Center 
University of Arizona 
Tucson Arizona 85721 


September 1987 


• (NASA-CB-181236) IKE AESCLUlf E ADICMETBIC N87-27969 

CALIE RATION OE IKE ADVA EC EE VIEY KIGfi 
EESCLCTICN BAEICBEIEE Semiannual Status 

Eeport (Arizona Gniv.) 27 p Avail: Nils Unclas 

EC AC3/BE AC 1 CSCL 14B G3/35 0097165 



1. Introduction 


An increasing number of remote sensing investigations require radiometrically 
calibrated imagery from the NOAA AVHRR sensors. Although a pre-launch 
calibration was done, there is no capability for monitoring any changes in calibration 
in-orbit for the visible and near infrared spectral bands. Hence, the possibility of 
using the reflectance-based method developed at White Sands for in-orbit calibration 
of Landsat TM data (Slater et al., 1987) to calibrate the AVHRR sensor has been 
under investigation. 

The purpose of this document is to provide a status report on research activities 
under NASA grant NAG 5-859. Section 2 describes the methodology and presents 
early results obtained for the NOAA -9 AVHRR. Section 3 summarizes the current 
status and plans for future activities. 

Those involved in the preparation of this report were P. M. Teillet 1 , Yalan Mao 
and P. N. Slater. Yi Ding helped with the data reduction. Others involved in data 
collection at White Sands and Edwards Air Force Base, partially in support of this 
activity, are R. J. Bartell, S. F. Biggar and R. P. Santer from the Optical Sciences 
Center and R. D. Jackson and M. S. Moran from the USDA Water Conservation 
Laboratory in Phoenix. 

2.0 Methodology and Early Results for NOAA-9 AVHRR 

The reflectance-based approach to radiometric calibration at White Sands is 
strongly dependent on the accurate characterization of ground reflectance (among 
other things) at the time of satellite overpass. Suitable reflectance-factor 
measurements have been made over terrain areas corresponding to numerous Landsat 
TM or SPOT HRV pixels (Slater et al., 1987; Jackson et al., 1987). However, this 
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approach becomes impractical for the calibration of NOAA AVHRR image data 
with pixel dimensions of 1.1 km by 1.1 km or greater. One alternative is to acquire 
AVHRR and TM imagery on the same day and take advantage of the detailed TM 
calibration results to effect a calibration of the AVHRR sensor. 

On August 28, 1985, ground reflectance-factor measurements in TM bands, and 
atmospheric data were collected at Chuck Site, White Sands, New Mexico. The 
methodology and results for reflectance-based calibrations of the Landsat-5 TM on 
this date and four others have been presented elsewhere (Slater et al., 1986). These 
show a ± 2.8% uncertainty in precision (la RMS standard deviation) for the six 
solar-reflective bands over the five dates. Table 1 lists the results for August 28, 
1985. The TM band 4 values differ slightly from those reported previously (Slater et 
al., 1986) because of a refinement in the way water- vapor absorption is handled. 

Given the calibration of TM bands 3 and 4 and NOAA-9 AVHRR imagery in 
bands 1 and 2 on the same day, the prediction of radiance at the AVHRR sensor 
may be divided into four steps: 

(i) determination of ground reflectance (using TM bands 3 and 4) averaged over 
as large an AVHRR pixel area as possible (essentially 2 pixels by 2 lines) in the 
alkali flat region of White Sands; 

(ii) spectral redistribution of ground reflectance into seven narrower bands 
encompassing AVHRR bands 1 and 2; 

(iii) prediction of radiance at the AVHRR sensor in the narrower bands using 
measured atmospheric parameters and radiative transfer computations; 

(iv) reconstitution of predicted radiances into spectral radiance values in AVHRR 
bands 1 and 2, and expression of the absolute radiometric calibration in terms of 
counts per unit radiance. 


Table 1. Calibration of the Thematic Mapper at White Santis on August 28, 1985. 
(*With additional gaseous transmittance due to H 2 O, O 2 , CO 2 ) 


Solar zenith angle Z: 35.954 

Solar distance in AU: 1.0098 

Junge size distribution: 3.77 

Aerosol size range: 0.02 to 5.02 n m 
Refractive index: 1.54 - 0.0 li 

Time of overpass: 10:08.4 MST 
Calculated visibility 200 km 


Latitude: 32 deg 55 min 
Longitude: 106 deg 22 min 
Elevation: 1196 m 

Pressure: 658.1 mm Hg or 877.42 mbars 
Temperature: 29.5 deg C 
Relative humidity: 21% 

Nadir viewing angle: 5 deg. 


Thematic Mapper bands 

3 

4 

4* 

Central wavelength pm 

0.6607 

0.8382 

0.8382 

Tau Mie 

0.0588 

0.0386 

0.0386 

Tau Rayleigh 

0.0403 

0.0155 

0.0155 

Tau ozone 

0.0141 

0.0016 

0.0016 

Tau water vapor 

0.0000 

0.0341 

0.0341 

Tau carbon dioxide 

0.0000 

0.0000 

0.0000 

Single scattering albedo 

0.8133 

0.7824 

0.7824 

Spectral reflectance 

0.6018 

0.6422 

0.6422 

Eo across band in W/m2.jun 

1545.0 

1042.8 

1042.8 

Average image digital counts 

203.31 

166.14 

166.14 

Preflight cal gains 

10.203 

10.821 

10.821 

Preflight cal offsets 

1.8850 

2.2373 

2.2373 

IC cal gains for 28 Aug 1985 

9.25 

10.299 

10.299 

IC cal offsets for 28 Aug 1985 

3.1170 

2.7320 

2.7320 

Normalized code radiance 

0.1471 

0.1496 

0.1538 

Code TM L in W/m2.sr.^m 

222.84 

153.02 

157.36 

Spectral L from preflight cal 

197.42 

151.47 

151.47 

Spectral L from IC cal 

216.42 

158.66 

158.66 

% (Code-Pre)/Pre 

12.9 

1.0 

3.9 

% (Code-IC)/IC 

3.0 

-3.6 

-0.8 

Counts per unit radiance 

0.912 

1.086 

1.056 

***************************************************************************** 

Case for Rayleigh Atmosphere 

Normalized code radiance 

0.1555 

0.1655 


Code TM L in W/m2.sr./im 

235.57 

169.28 


Counts per unit radiance 

0.863 

0.981 



***************************************************************************** 


Case for No Atmosphere 


Normalized radiance 
TM L in W/m2.sr.^m 
Counts per unit radiance 


0.1551 0.1655 

234.90 169.20 

0.866 0.982 


***************************************************************************** 



2.1 Determination of Ground Reflectance from TM Imagery 


Since the TM and AVHRR sensors acquired images with very different obser- 
vation geometries (TM with near-nadir viewing in a descending orbit, AVHRR with 
significant off-nadir viewing in an ascending orbit), a geometric registration 
procedure was used to match the two images. From the superimposed images 
(Figures 1-4), a relatively uniform area of two by two AVHRR pixels was selected 
in the alkali-flat region of White Sands. The digital counts for this area and for the 
corresponding area in the TM imagery were extracted and averaged. The TM values 
are presented in Table 2, together with radiance values based on the White Sands 
calibration from Chuck Site (Code) and based on the TM Internal Calibrator (IC). 
Results are also given for an area s imil ar in size in the dunes area at White Sands. 

In order to relate the TM radiance values obtained for the AVHRR test area to 
ground reflectance factors, a series of atmospheric computations were carried out 
using the Herman radiative transfer code (Herman and Browning, 1975) and the 
atmospheric input parameters listed in Table 1. With ground reflectance values 
ranging from 0.00 to 0.99 at the input, corresponding radiance values at satellite 
altitude were generated by the code. For both TM bands 3 and 4, the resulting 
relationship between predicted radiance and ground reflectance is quite linear. These 
linear relationships were then used to obtain reflectance factors from the code 
radiance values derived from TM bands 3 and 4 for the AVHRR test area (Table 2). 

Table 2. Ground reflectances for the AVHRR test area from TM imagery. 

Average L(IC) L(CODE) Central 

TM Band Area DC (W/M2.SR./JM) (W/M2.SR./iM) Reflectance Wavelength 


3 

Alkali Flat 

172.3 

182.9 

188.3 

0.510 

0.6607 

3 

Dunes 

239.9 

256.0 

263.6 

0,720 

0.6607 

i 

Alkali Flat 

144.2 

137.3 

132.4 

0.555 

0.8382 

4 

Dunes 

195.6 

187.2 

180.5 

0.760 

0.8382 
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Figure 1. Digital counts of NOAA-9 AVHRR channel 1 at White Sands on August 28, 
1985. The higher values in the upper portion correspond to the dunes 
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Figure 2. Digital Counts of NOAA-9 AVHRR channel 2 at White Sands on August 28, 
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Figure 3. 


Digital counts of Landsat-5 TM channel 3 for alkali-flat region of 
White Sands on August 28, 1985. 

144 U4 147 172 173 174 170 171 177 173 173 171 174 173 170 170 

141 U3 147 172 171 ITS 179 177 171 171 173 173 174 172 172 144 

142 US 144 140 171 171 178 177 177 173 173 174 144 144 144 141 

144 US 147 147 148 173 173 173 US 173 174 174 173 148 138 141 

141 171 172 171 172 173 174 171 144 172 174 174 173 173 143 142 

14» 171 174 172 174 177 177 174 173 144 144 U7 173 172 144 149 

171 171 171 172 172 174 177 174 173 173 174 Uf 144 172 144 142 

171 ITS 171 172 171 173 ITS 174 174 174 174 171 141 170 144 140 

141 170 173 170 177 177 174 174 174 174 173 ITS 174 173 174 173 

144 US 174 173 102 103 ITS 173 174 173 173 174 173 ITS 173 173 

170 171 174 ITS 173 170 170 174 173 ITS 177 172 173 171 173 173 

171 172 174 ITS 173 ITS 113 113 102 102 174 172 170 144 170 177 

172 172 172 173 172 174 104 104 103 102 174 170 ITS 173 U3 171 

170 173 143 143 172 175 182 187 180 148 143 170 145 143 174 172 

150 143 143 145 143 174 173 184 182 US 170 170 14? 143 173 172 


154 154 147 171 170 174 173 184 1S4 130 174 143 148 173 173 143 
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Figure 4. Digital counts of LANDSAT-5 TM channel 4 for alkali-flat reigion of 
White Sands on August 28, 1985. 
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, ; 
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2.2 Spectral Redistribution of Ground Reflectance 


Since NOAA-9 AVHRR spectral bands 1 (0.570-0.700 /an) and 2 (0.714-0.983 
/im) are not narrow bands (Figure 5) and the Herman radiative transfer code assumes 
monochromaticity, they were subdivided into small segments. Channel 1 was divided 
into three bands centered at 595, 635, and 680 nm, whereas channel 2 was divided 
into four bands centered at 760, 840, 900, and 960 nm. Before running the Herman 
code at these wavelengths, the reflectance values obtained in the previous section 
(2.1) were adjusted in two ways. 

Firstly, at the NOAA-9 satellite overpass time of 21:27 Universal Time, the 
solar zenith angle was 39.85 degrees, whereas at the Landsat 5 satellite overpass time 
of 17:08 Universal Time, the solar zenith angle was 35.95 degrees. Moreover, the 
viewing nadir angle was 23.62 degrees for the AVHRR sensor and 5 degrees for the 
TM sensor. Thus, corrections for non-lambertian effects were applied to the TM 
band 3 and 4 reflectance factors on the basis of bi-directional reflectance 
measurements made for the gypsum surface at White Sands on March 15, 1986 
(Figure 6). 

Secondly, the reflectance factors derived in the previous section were for TM 
bands 3 and 4, but ground reflectances for input to the Herman code were required 
at seven other wavelengths. Fortunately, off-nadir ground reflectance measurements 
with a Barnes Modular Multispectral 8-channel Radiometer 2 (Robinson et al., 1979) 
are available for approximately the same solar zenith as that for the AVHRR 
overpass on August 28, 1985. Those measurements were used by Ray Jackson to 
derive an equation relating the surface reflectance at Chuck Site to wavelength in 


2 Trade names and company names are included for the convenience of the reader 
and imply no endorsement of the product or company by the University of Arizona 
or the United States Department of Agriculture. 
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0 micrometers (Table 3): 

p(A) * a 0 + a x A + a 2 A 2 + a 3 A 2 + a 4 A 4 . 

With this equation and coefficients from Table 3 for the case of 24.6 degrees nadir 
9 view angle, reflectances were calculated at the seven wavelengths of interest (Table 

4). These values were also adjusted to take into account the difference in surface 
reflectance between Chuck Site and the alkali-flat region, the latter being darker in 
# the spectral domain of interest. More specifically, the three reflectances in the 

AVHRR band-1 region were shifted downward by the difference between p( A = 
0.6607) from Jackson's equation and the relevant reflectance value (0.510) from 
® Table 2. Similarly, the four reflectances in AVHRR band-2 region were shifted 

downward by the difference between p ( A * 0.S382) and the relevant reflectance value 
(0.555) from Table 2. 


Table 3. Surface Reflectance as a Function of Wavelength for Different View 
Angles. 


Date 

Solar Zenith 
Angle 
(Degree*) 

View Nadir 
Angle 
(Degrees) 

*0 

a i 

a 2 

a 3 

a 4 

08 Mar 86 

40.10 

31.2 

-0.84395 

6.00165 

-9.18091 

6.2899 

-1.6215 

15 Mar 86 

41.01 

24.6 

-1.85425 

10.69285 

-18.28977 

14.1966 

-4.1502 

20 Mar 86 

38.62 

18.4 

-1.64884 

9.53550 

-16.16274 

12.5705 

-3.7022 


Table 4. Reflectances for the Seven Narrower Bands Subdividing the AVHRR 
Bands. 


Wavelength A 
(moi) 

Reflectance p( A) 
(Chuck Site) 

Shift 

(See Text) 

Reflectance 
(Alkali Flat) 

0.595 

0.5032 

0.0102 

0.4930 

0.635 

0.5210 

0.0102 

0.5108 

0.680 

0.5362 

0.0102 

0.5260 

0.760 

0.5555 

0.0193 

0.5362 

0.840 

0.5706 

0.0193 

0.5513 

0.900 

0.5810 

0.0193 

0.5617 

0.960 

0.5903 

0.0193 

0.5710 
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Figure 5a. Relative Spectral Response for NOAA-9 AVHRR Channel 1 (Kidwell 
1986). 
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Figure 5b. Relative Spectral Response for NOAA-9 AVHRR Channel 1 (Kidwell 
1986). 



2.3 Prediction of Radiance at the AVHRR Sensor 


Atmospheric parameters and ground reflectances for the seven wavelengths were 
input to the Herman radiative transfer code. Output from the code includes 
radiance at the entrance aperture of the AVHRR sensor, normalized for unity exo- 
atmospheric solar irradiance. For use in the code, the atmosphere is divided into a 
sufficient number of plane-parallel layers such that changes within each layer are 
due only to single-scattering processes. The Gauss-Seidei iterative technique is used 
to solve the equation of radiative transfer. Upon convergence, all multiple scattering 
effects have been taken into account. Values of 5.02, 0.02, and 0.04 /xm were used 
for the maximum and minimum radii and incremental step size, respectively, for the 
aerosols. A vertical aerosol distribution as measured by Elterman (1966) was 
assumed. The aerosols were given a refractive index of 1. 54-0.0 li, an average value 
for the region as measured by Jennings et al. (1978). 

Other important inputs to the Herman code include Rayleigh, aerosol, and ozone 
optical depth values. The Rayleigh optical depth was determined from a knowledge 
of the barometric pressure and wavelength. The total optical depths were 
determined from the slopes of Langley plots in which the log voltages from solar 
radiometers were plotted against air masses. In spectral regions unaffected by 
absorption, the Mie optical depth at any wavelength was determined by subtracting 
the Rayleigh from the total optical depth at that wavelength. An optical depth 
versus wavelength curve was fitted through the points that spanned the absorption 
region due to ozone (Figure 7). The differences between the values on this curve at 
a given wavelength and the total minus Rayleigh value at the same wavelength gave 
the ozone optical depth at that wavelength. A Junge radial size distribution (Junge, 
1963) was v assumed for the aerosols. The Junge v parameter is defined in the 
equation 
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Figure 6. Reflectance factor variations as a function of view angle for two 
solar zenith angles on March 15, 1986. The curves are for TM bands 

1 to 4 with the bands in order starting with band 1 at the bottom in 
each case. 
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Figure 7. Optical depth versus wavelength plot used to determine ozone optical 
depth. 

Aug. 28. 1985 
Chuck Si to 




where N is the number of particles, r is their radius, and C is a constant. The value 
of v , needed to find the aerosol phase function, was determined from the slope of 
the log versus log A curve. 

The results are summarized in Tables 5 and 6. The code spectral radiance is the 
product of the exo-atmospheric irradiance and the normalized radiance determined 
from the Herman code. The exo-atmospheric solar irradiance data are those 
recommended by Frohlich and published by Iqbal (1983). They represent a carefully 
edited combination of results published by Neckel and Labs, Thekaekara, Arvesen, 
and others. The values were adjusted to yield an integrated value of 1367 Wm -2 , 
the solar constant as proposed by the World Radiation Center. 

2.4 Reconstitution of Predicted Radiances into the AVHRR Bands 

In order to combine the code spectral radiance results for the narrower bands, 
weighting coefficients were determined from the areas under selected intervals of the 
relative spectral response profiles for each AVHRR band (Figure 5). For example, 
the weighting coefficient for 595 nm is given by 

-610 

R(A)dA 

J530 

Tsoo ’ 

R(A)dA 

J530 

The code spectral radiance in each AVHRR band is then the weighted sum of code 
spectral radiances for the narrower bands (Tables 5 and 6). The central wavelengths 
listed for the AVHRR bands were calculated by the moments method (Palmer, 


Table 5. Calibration of NOAA-9 AVHRR Band 1 at White Sands on August 28, 1985. 


Solar Zenith Angle Z: 39.85 

Solar Distance in AU: 1.0098 

Junge Size Distribution: 3.77 

Aerosol Size Range: 0.82 to 5.02 pm 

Refractive Index: 1.54 - 0.0 li 

Time of Overpass: 21:27 (universal time) 

Calculated Visibility: 200 km 


Location: Alkali Flat 
Latitude: 32 degree 47 min. 

Longitude: 106 degree 22 min. 

Satellite Zenith Angle: 23.62 Degree 
Satellite Azimuth Angle: 255.00 Degree 


Wavelength in nm 
Bandwidth 

Tau Mie 
Tau Rayleigh 
Tau Ozone 
Spectral Reflectance 
Exoatmospheric E0 (W/m2./xm) 
Normalized Code L 
Spectral Radiance 
Weighting Coefficients 
Weighted Spectral Radiance 


Total Gas Transmittance: 

Ozone Transmittance: 

Total Gas Transmittance Except Ozone: 
Radiance Predicted at the Sensor (Pcode): 
%(Pcode-pre)/pre: 

%(5S-Pcode)/Pcode: 

Counts Per Unit Radiance: 


595 

635 

680 

530 - 610 

610-660 

660-800 

0.0711 

0.0633 

0.0561 

0.0618 

0.0474 

0.0359 

0.0338 

0.0219 

0.0102 

0.4930 

0.5110 

0.5260 

1785 

1630 

1472.5 

0.1079 

0.1151 

0.1220 

188.86 

184.10 

176.20 

0.2912 

0.3677 

0.3411 

54.988 

67.686 

61.256 


0.6400 

183.94 

367.25 (10-bit scale) 

172.25 
6.78 
1.9967 
2.1321 
0.5100 
0.4730 

178.47 

3.61 

0.9400 

0.9550 

0.9843 

181.04 

5.10 

-1.42 

2.0285 


Central Wavelength for (/un): 

Code Spectral Radiance (W/m2./im.sr): 
Average Image Digital Counts (2 by 2 pixels): 
Spectral Radiance from Preflight 
Calibration (W/m2./im.sr): 

%(Code-Pre)/Pre: 

Code Counts Per Unit Radiance: 

Preflight Counts Per Unit Radiance: 

Average Reflectance from TM3 (8/28/85): 
Apparent Reflectance R (5S): 

Spectral Radiance from 5S 
(L=R*cosZ* 1 544/3. 14159 W/m2.Mm.sr): 
%(5S-Pre)/Pre (Radiance): 


Table 6. Calibration of NOAA-9 AVHRR Band 2 at White Sands on August 28, 1985 


Solar Zenith Angle: 39.85 

Solar Distance in AU: 1.0098 

Junge Size Distribution: 3.77 

Aerosol Size Range: 0.02 to 5.02 pm 

Refractive Index: 1.54 - 0.0 li 

Time of Overpass: 21:27 (universal time) 

Calculated Visibility: 200 km 


Wavelength in nm 760 

Bandwidth 680-800 

Tau Mie 0.0459 

Tau Rayleigh 0.0229 

Tau Ozone 0.0040 

Spectral Reflectance 0.5360 

Exoatmospheric E0 (W/m2.pm) 1222.5 

Normalized Code L 0.1265 

Spectral Radiance 151.71 

Weighting Coefficients 0.3403 

Weighted Spectral Radiance 51.624 


Central Wavelength (pm): 

Code Spectral Radiance (W/m2.pm.sr): 
Average Image Digital Counts (2 by 2 pixels): 
Spectral Radiance from Preflight Calibration 
(W/m2.pm.sr): 

%(Code-Pre)/Pre: 

Code Counts Per Unit Radiance: 

Preflight Counts Per Unit Radiance: 

Average Reflectance from TM4 (8/28/85): 
Apparent Reflectance R (5S): 

Spectral Radiance from 5S 
(L=R*cosZ*999 .2/3. 14159 W/m2.pm.sr): 
%(5S-Pre)/Pre (Radiance): 

Total Gas Transmittance: 

Ozone Transmittance: 

Total Gas Transmittance Except Ozone: 
Radiance Predicted at the Sensor (Pcode): 
%(Pcode-pre)/pre: 

%(5S-Pcode)/Pcode: 

Counts Per Unit Radiance: 


Location: Alkali Flat 
Latitude: 32 degree 37 min 
Longitude: 106 degree 22 min 
Satellite Zenith Angle: 23.62 Degree 
Satellite Azimuth Angle: 255.00 Degree 


840 

900 

960 

800-870 

870-950 

950-1160 

0.0385 

0.0340 

0.0304 

0.0153 

0.0116 

0.0089 

0.0016 

0.0000 

0.0000 

0.5510 

0.5628 

0.5710 

1020.0 

913.0 

771.0 

0.1311 

0.1344 

0.1368 

131.15 

120.37 

103.41 

0.2601 

0.2728 

0.1268 

34.118 

32.833 

13.115 


0.8460 

131.69 

327.31 (10-bit scale) 

96.633 

36.28 

2.4855 

3.3872 

0.5550 

0.4660 

113.79 

17.75 

0.8490 

0.9980 

0.8507 

112.03 

15.93 

1.57 

2.9217 


A final adjustment was applied to the code spectral radiances to correct for 
gaseous absorption. The Herman code takes into account ozone optical depth but 
does not consider H 2 0, C0 2 , and 0 2 transmittance. In AYHRR band 2, there is 
significant absorption around 0.94 /mi due to water-vapor and around 0.76 /xm due 
to oxygen. Thus, the "5-S" atmospheric program of Tanre et al. (1985) was run to 
obtain the total gaseous transmittance due to the four gases (H 2 0, 0 3 , C0 2 , 0 2 ). 
The final radiance values predicted at the AYHRR sensor are the product of total 
gas transmittance except ozone and the code spectral radiance (Tables 5 and 6). 

2.5 Discussion of Calibration Results 

Compared to the pre-launch calibration (from February, 1980) for the NO A A -9 
AVHRR sensor (Smith, 1987), the results indicate response degradations of 5.1% and 
15.9% in AVHRR bands 1 and 2, respectively. The outcome of this analysis differs 
significantly with the results of a high-altitude aircraft calibration of the NOAA-9 
AVHRR carried out over White Sands on August 26, 1985 (Table 7). According to 
Smith et al. (1987a), the aircraft calibration results show no degradation in channel 1 
and 6.3% degradation in channel 2. A more detailed analysis of the same data set 
(Smith et al., 1987b) indicates degradations of +2 ± 5% for channel 1 and -2 ± 5% 
for channel 2, i.e. no significant change. Flight data from October and November 
1986 have also been analyzed (Smith et al., 1987a), indicating decreases in sensor 
response of 14.3% and 23.2% in AVHRR bands 1 and 2, respectively. 

The work of Frouin (1987) and Frouin and Gauthier (1987) on AVHRR 
calibration should also be mentioned. Using standard parameters, the dunes area at 
White Sands and dark space digital counts, but without reference to any near- 
coincident ground or atmospheric measurements, they monitored the gain of several 
satellite sensors as a function of time. For the NOAA-9 AVHRR in the same 
October- November 1986 time frame studied by Smith et al. (1987a), the results from 
this completely different approach are quite close, with degradations of 12.9 ± 8% 



for channel 1 and 23.7 ± 8% for channel 2 (Frouin, 1987). Frouin and Gauthier 
(1987) monitored the NOAA-7 AVHRR sensor from September 1983 to December 
1984 using the same technique, encompassing the analysis of at least 100 pixels in 
each of twelve AVHRR scenes during that time period. Their results show no 
systematic changes in gain during that time, but the response degradation with 
respect to the pre-launch calibration (Kidwell, 1986) was 14 ± 8% for AVHRR band 
1 and 15 ± 8% for band 2 (Table 7). 


Table 7. Change in sensor response with time since launch (in percent), with 
respect to pre-launch calibration coefficients from Smith (1987) for 
NOAA-9, and from Kidwell (1986) for NOAA-7. 


Months 

AVHRR 






Since 

Channel 

NOAA-9 

NOAA-9 

NOAA-9 

A A A 

NOAA-7 

Launch 

Number 

AVHRR 1 

AVHRR 2 

AVHRR 3 

AVHRR 4 

AVHRR 5 


1 

- 5.1% 

0% 

-2 ±5% 



8 

2 

-15.9% 

-6.3% 

+2 ± 5% 




1 


-14.3% 


-12.9 ± 8% 


22 

2 


-23.2% 


-23.7 ± 8% 



1 





-14 ± 8% 

28 

2 





-15 £ 8% 


1 -14 ± 8% 


* Thi3 Study 

2 Smith et al., 1987a 

3 Smith et al., 1987b 

4 Frouin, 1987 

3 Frouin & Gauthier, 1987 




It is worth noting that the change between August 1985 and October/November 
1986, according to Smith et al. (1987a), is considerably larger than that between 
launch (December, 1984) and August 1985. This trend differs from what has been 
observed with other sensors such as CZCS and TM, where the degradation was 
greatest in the months just after launch and the rate of change decreased thereafter. 

The reasons for the discrepancy between our ground-based reflectance approach 
and NOAA's aircraft-based radiance approach for August 1985 are not known, 
although the two methods are completely different in almost all respects. It is hoped 
that the analysis of additional data sets by both groups will clarify the situation. 

As a final comment in this section, Frouin (1987) has pointed out that there are 
several different values for the NOAA-9 AVHRR pre-launch gain coefficients 
available (Table 8). Thus, when monitoring and reporting changes in sensor 
responsivity with time, it is important to base results on a consistent set of pre- 
launch values. 


Table 8. Various Versions of NOAA-9 AVHRR Pre-launch Gain Coefficients 
Expressed in Radiance (W/m^sr/im) per unit count. 


Source 

NOAA-9 
AVHRR Band 1 

NOAA-9 
AVHRR Band 2 

NASA/GSFC 

0.5243 

0.3286 

Smith (1987) 

0.5246 

0.3363 

Kidwell (1986) 

0.5191 

0.3724 

3.0 Current Status and Plans 


The current status of the AVHRR calibration studies, including future plans, 
can be summarized as follows. 


(1) A first pass through the analysis of the data set for August 28, 1985 has been 
completed. Preliminary results indicate significant degradation in the response of 


both AVHRR bands 1 and 2. These results are currently in the process of being 
verified. 

(2) Additional data sets involving AVHRR data, same-day imagery from a higher- 
resolution sensor, and ground measurements are being prepared (Table 9). In 
particular, the next data set to be analyzed consists of Landsat-5 TM and NOAA-9 
AVHRR image data acquired on March 27, 1987 at White Sands together with 
surface reflectance and atmospheric data collected in the usual manner for Chuck 
Site. In general, it has not proven to be easy to obtain AVHRR and TM data sets 
on the same field sortie day at White Sands, with the AVHRR nadir view angle less 
than about 40° and the weather remaining clear for both overpasses. Another 
problem which was not anticipated is the difficulty in obtaining archival AVHRR 
data. A complete archive does not exist at NASA, NOAA, Scripps or the EROS 
Data Center, and this has severely limited out ability to perform retrospective 
comparisons of AVHRR calibrations. We have also encountered problems with the 
reliability of data acquisitions from both TM and HRV. In the TM case it appeared 
to be due to a breakdown in communications between EQSAT and the Prince Albert 
receiving station and in the HRV case it was due to an error in the input of the 
sensor pointing angle. In both cases it involved the loss of potentially valuable data 
together with the considerable waste of time and effort of field teams numbering 
about ten people in each case. 

(3) During the next six months, field sorties are being planned to coincide with 
TM and/or HRV overpasses and with same-day AVHRR overpasses for which nadir 
view angles do not exceed 40 degrees. It is hoped that two or three additional study 
cases can be acquired and analyzed for the purpose of AVHRR radiometric 
calibration, including the NOAA- 10 AVHRR instrument. 


Table 9. Status of Sensor Data Reduction (August 1987) 

Overpass 

Date Place Sensor Time (UT) SZA SAA VAG AAG Comments 


28 Aug 85 

WSMR 

TM (Landsat 5) 
AVHRR (NOAA 9) 

17:08:24 

21:27:00 

14 Oct 86 

EAFB 

AVHRR (NOAA 9) 
AVHRR (NOAA 10) 
AIS 


27 Mar 87 

WSMR 

AVHRR (NOAA 10) 
TM (Landsat 5) 

15:15:45 

17:01:03 

28 Mar 87 

WSMR 

AVHRR (NOAA 10) 
HRV-1 (SPOT 1) 

12:54:05 

18:18:21 

04 May 87 

EAFB 

AVHRR (NOAA 9) 
AVIRIS 

22:29:54 

05 May 87 

EAFB 

AVHRR (NOAA 9) 
AVIRIS 

AVHRR (NOAA 10) 

22:19:03 

16:10:14 

17 Jul 87 

WSMR 

AVHRR (NOAA 10) 

14:46:32 



HRV-? (SPOT 1) 

17:39:06 


36 

124 

5.0 


Preliminary 

40 

242 

23.6 

255.0 

Analysis Completed. 


62 

107 

21. SW 

284.6 

CCT from NOAA read. 

43 

129 

0 

0 

CCT from EOSAT read. 

65 

104 

12.5E 

31.2W 

102.4 

278.7 

CCT from NOAA read. 

DC values received from 
CNES. 

41 

252 

15.3E 

81.7 

CCT from Scripps read. 

39 

250 

31.3E 

81.7 

CCT from NOAA read. 
CCTs received from JPL. 

52 

97 

28.2W 

279.4 

CCTs unavailable. 

59 

183 

32.3E 

99.2 

CCT requested from EDC 
and Scripps. 

23 

113 

30.2E 

99.0 

DC values requested from 
CNES. 


Latitude 

WSMR = White Sands Missile Range 32° 55’ 
EAFB = Edwards Air Force Base 35° 00’ 


Longitude 
106° 22’ 
117° SO’ 


Elevation 
1196 m. 
694 m. 


SZA = solar zenith angle 

SAA = solar azimuth angle 

VAG = viewing angle from ground to sensor 

AAG = azimuth angle at ground 
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(4) Other measurement and analysis possibilities under consideration include the 
collection of ground-based data at Chuck Site at the time of the AVHRR overpass, 
and the use of the extended black lava area at White Sands in both TM and AVHRR 
imagery as a dark target for estimating atmospheric parameters. 

(5) The current methodology makes use of the Herman code, supplemented by the 
"5-S" code and actual ground-based measurements of atmospheric conditions and 
surface reflectance, to predict spectral radiance values at the AVHRR sensor. The 
approach takes advantage of a reasonably precise calibration of TM image data 
acquired on the same day to derive an absolute calibration for the AVHRR imagery. 
It is of interest to learn to what extent the calibration process can be simplified, 
possibly even to the point where no ground measurements are required, for relative 
intercalibration purposes. Initially, this investigation will involve the "5-S" 
atmospheric code modified to give results in an inverse mode, accepting apparent 
reflectance at satellite altitude and computing the ground reflectance. In this way, 
TM imagery, assumed to be well calibrated independently, can be used with "5-S" to 
generate ground reflectances for the alkali-flat region of White Sands. With suitable 
adjustments for non-lambertian effects and wavelength dependence, these ground 
reflectances can then be used with "5-S" in the forward mode to predict radiances 
at the AVHRR sensor on the same day. It may be that the unknown atmospheric 
conditions will not affect the calibration significantly because of the two-way use of 
the "5-S* code. More specifically, the sensitivity of this procedure to uncertainties in 
key atmospheric parameters will be examined. 

If results from this approach compare favorably with the more detailed 
methodology and are not overly sensitive to atmospheric conditions, the implication 
is that a reasonable calibration of satellite sensors with no on- board calibration 
capability may be possible by transfer without the necessity of making ground-based 
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measurements. In this way, it would be relatively straightforward to monitor 
occasionally the status of the AVHRR sensor's radiometric response. 

4.0 Concluding Remarks 

The early results of an absolute radiometric calibration of the NOAA-9 AVHRR 
sensor indicate significant degradations in the response of bands 1 and 2 compared to 
pre-launch values. The results are currently in the process of being verified and it 
may be that refinements of the methodology will be in order as additional data sets 
are analyzed. Consequently, the reported results should not be considered conclusive 
for the time being. Nevertheless, the Landsat TM calibration used in this approach 
is known to be very precise and the Herman radiative transfer code, supplemented 

V*V ttiA Tnr o«*CAfYn« ie voltohla wall XVia avtont +r\ 

other steps in the analysis procedure give rise to uncertainties in the results is 
currently under investigation. Particular attention is being given to the geometric 
matching of the AVHRR and TM imagery, as well as to the spectral redistribution 
procedure. 

By taking advantage of a reasonably precise calibration of TM imagery acquired 
on the same day as the AVHRR data at White Sands, a promising approach to the 
in-orbit calibration of AVHRR sensors is being developed. Current efforts involve 
primarily the examination of additional test cases and the investigation of possible 
simplifications in the procedure through judicious use of atmospheric models. 
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